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Using the Thomas-Fermi approximation, we show that an 
interacting two dimensional electron gas may be described in 
terms of fractional exclusion statistics at zero and finite tem- 
peratures when the interaction has a short-range component. 
We argue that a likely physical situation for this phenomenon 
to occur may exist in two dimensional quantum dots. 
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Fractional exclusion statistics or the generalised exclu- 
sion principle was first proposed by Haldane in the 
context of excitations in spin chains. Experimentally, the 
best evidence comes from recent neutron inelastic scat- 
tering experiment |^ on the compound KCuF^, which 
is a one-dimensional Heisenberg antiferromagnet above 
AO'^K. The observed inelastic scattering is best fitted by 
spinon excitations in a spin chain whose pairwise interac- 
tion falls off as the inverse square of the lattice distance 
Q. The dynamic correlation function for such a system 
has been calculated by Haldane and Zirnbauer The 
concept of fractional exclusion statistics has been gen- 
eralised to the case of a gas of particles Q defined by a 
distribution function 1^,^ that allows for partial or multi- 
ple occupancy of a single-particle state. In principle, the 
statistics is applicable to particles in any spatial dimen- 
sion, but most known examples are mathematical mod- 
els in one dimension with pairwise inverse-square 
interaction. The first calculation for a two-dimensional 
realistic system in this context was done by Johnson and 
Canright who demonstrated, by exact diagonalisa- 
tion of a small number of interacting electrons, that the 
bulk excitations in FQHE liquids exhibit Haldane statis- 
tics. In this paper, we show that under certain condi- 
tions, a two-dimensional interacting electron gas in its 
ground state may exhibit this statistics. The conditions 
are shown to be favourable for electrons in a quantum 
dot. In this case, it is shown that the dominant effect 
of the interaction may be incorporated in the fractional 
statistics of the gas. If the residual interactions are ne- 
glected, then that the system also obeys Haldane statis- 
tics at finite temperature. This opens up the exciting 
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possibility that the bulk properties of a mesoscopic two- 
dimensional system may be understood by regarding it 
as an almost ideal fractional statistics gas confined in a 
potential well. 

The claims made in this paper are based on the 
Thomas- Fermi (TF) method [Q. Being a mean-field 
method, it cannot reproduce two-body correlations, but 
is successful in giving a good estimate of bulk properties 
like the ground-state energy and the single-particle spa- 
tial density. It has previously been applied with success 
to atoms llj] , nuclei ||l^ , and metal clusters [|l6| . In two- 
dimensions, TF yields an accurate approximation to the 
total energy of a many-anyon system |17|. For an ideal 
gas obeying the generalised exclusion statistics, TF cal- 
culation has been shown to yield the exact answer for the 
energy in the large-N limit ||l^ . It is therefore reasonable 
to expect that the method gives meaningful answers. We 
start by constructing the energy density functional for 
the ground state energy of a system of interacting spin- 
half fermions. Consider the N-fcrmion Hamiltonian in 
two dimensions: 

^ N N 

i—1 i—1 j<k 

where Vi is a one-body confining potential whose specific 
form is not crucial at present and V2 is the two body 
potential which is repulsive. In a mean-field theory, the 
expression for the energy at zero-temperature is given by. 
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2m 
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where p{r) is the spatial single-particle density, T(r) is the 
kinetic energy density and pir^ r') is the density matrix. 
In the above we have taken into account the effect of both 
direct and exchange terms in the interaction energy. The 
factor 1/2 is the correction due to over-counting of pairs. 
The constant C is determined by the spin polarisation of 
the gas : for unpolarised electrons, it is 1/2, whereas for 
a fully polarised system, it is 1. For arbitrary polariza- 
tion P = '^^'^^^ , where N± is the number of up or down 
spins, the factor C = 
malized such that TV 



2 . The spatial density is nor- 
J d?rp{r). In the Thomas-Fermi 
method, the kinetic energy density T{r) is itself expressed 



1 



in terms of the density p{r) and its gradients. The en- 
ergy and the density are determined self-consistently by 
a variational principle. In two-dimensions, the TF ex- 
pression for T(r) = 7rp^(r), taking into account the spin- 
degeneracy factor of 2. In this case, there is no gradient 
correction in the bulk up to 0{h'^). However, there are 
edge corrections when the sample is of finite size p9| . 

Next consider the energy due to the two body interac- 
tions. The matrix elements of the direct term is, 

< ij\V2\ij >= / p{ri)p{r2)V{\ri ~ f2\)d^rid^r2, 

(3) 

where the sum (here as well as in what follows) is over the 
occupied single-particle states only. The matrix elements 
of the exchange term is 

5]<*j1V^2|j«>- / \piri,r2)fVi\n-r2\)d\id\2, 

(4) 

where p{ri,r2) = Si V'j*(''i)'0i(''2)- At this stage, it is 
useful to perform the density-matrix expansion following 
Skyrme pO| . Defining r = ri— r2 and R = (ri-|-r2)/2 and 
expanding the density up to this order in r , we obtain 



pin) = piR + r/2) = p{R) + (r'.V)p + -(r.V)V + 



(5) 



The direct matrix element may then be written as, 
^ < ij\V2\ij > = J d\V2{r) J d^Rp\R) 

-\j d\ V2{t) J d^R{^p{R)f + . . . , (6) 

Similarly the density matrix p{'ri,r2) may be expanded 
up to second order in r about R, 

p{ri,r2) = ^**(^ + ^72)^. (i? - 72) 

i 
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and the exchange contribution to second order is given 
by, 

51 < *Jl^2|J^ j d\V2{r) J d^Rp'iR) 
_i J d^r V2ir) J d^RT{R)p{R) + . . . . (8) 



Here the kinetic energy density is defined as, 

i i 

(9) 

Often the kinetic energy density is defined either by the 
first term or by the second term in the above equation 
without the over all 1/2. What we naturally get in the 
expansion is an average of both these commonly used 
forms. We have computed each one of these forms ex- 
actly using harmonic oscillator wave functions for a few 
particles. While the the first and second terms show os- 
cillations around the smooth TF density, the definition 
given above almost precisely coincides with the TF den- 
sity even with as little as two particles. 

We note that the leading terms in both direct and ex- 
change terms are the same (proportional to p^). For spin- 
half fermions the interaction energy is given by, 



<^j\V2\3i>l (10) 



where is the spin projection. Summing over all par- 
ticle indices immediately gives a factor (1 -I- |P|)/2 for 
the exchange contribution, where P is the spin polariza- 
tion of the system. Therefore, if there is no other degree 
of freedom, or if the spins are all polarized, the contri- 
bution from the leading terms to the interaction energy 
vanishes as it happens in FQHE systems. However for 
the unpolarized 2-D electron systems there is a factor 
half for exchange contribution. Here we concentrate on 
the unpolarized case. Combining all the contributions 
the total energy of the system is given by. 



E 



fi 1 
— V(r) + Vi{r)p{r) + -p'{r)Mo 

+ \[{^p\r)-{yp{r)f]M2 + ..] , (11) 



where M„ — J (i^rV2(r)r" are the moments of the two 
body potential. Note that we obtain an expression simi- 
lar to the above if we use an expansion of the form [|l). 



(12) 



where b is the range of the potential and Cj are related to 
the j-th moment of the potential V2 as M2j — 2^^jlcjb^^. 

The spatial density is now determined by the variation 
S{E — p,N) — 0, where p is the chemical potential at 
zero temperature. The variation immediately gives the 
equation for the density. 
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In the large-N limit we expect the density in the bulk to 
be approximately constant. We can therefore neglect the 
derivative term in this limit. Further if the potential is 
extremely short-ranged, the term proportional to the sec- 
ond moment of the potential may also be neglected. (We 
will elaborate on these approximations shortly.) Then 
the density is given by, 



Poir) 



-0, 



(^-Vi(r)), r<ro 



r > ro, 



(14) 



where tq is the classical turning point defined by ii — 
Vi{ro) and 



a = 1 



m*Mo 



(15) 



is now the statistics parameter as we show below. In the 
effective range expansion (p2[), cq = Mq. The expression 



for pq in Eq.(14) may be interpreted as if the fermions in 
the one-body confining potential Vi are noninteracting, 
but that theyobey the generalised exclusion statistics for 
occupancy at zero temperature: 



a 

= 0, e > /i. 



(16) 



This may be easily seen as follows. For noninteracting 
fermions, the Thomas-Fermi density of states g{e) in an 
external potential Vi (r) is 



d^rd'^p 



5{e- 



P 

2m* 



-Vi{r) 



(17) 



The over-all factor of two on the right-hand side is due 
to the spin degeneracy. Using the new occupancies given 
by Eq.dl), we get 



N = 



1 /"^ 
a Jo 

1 f d'^rd'^p 

a J {2nhY 



2m' 



-Vi{r) 
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The function 9{y) — 1 for ?/ > 0, and zero otherwise. Now 
performing the p-integration immediately yields the total 
number of particles, with density po{r) given by Eg. (14). 
Indeed we have now the precise condition under which 
ideal exclusion statistics is realised within the framework 
of the Thomas-Fermi method. 

In the more realistic situation, the higher moments 
may not be neglected, and the system is a non-ideal frac- 
tional statistics gas. In the thermodynamic limit, we may 
write 



p(r) = po(r) 



^ 3m* M2 , , 

1 - -^72 — poin 

8n a 



(19) 



where po{r), given by Eq.(p^, is the density for the ideal 
FES case. Note that M2 = 4ci6^ where b is the range of 
the potential. The typical densities in two dimensional 
systems of interest is of the order of 10~^/A^. Using the 
values of m* = 0.067me, which is the effective electron 
mass in GaAs materials, and a > 1 (but not very large), 
it is easy to estimate that the second term becomes im- 
portant only for ranges of the order of lOOA or above. 
Another way to view the problem is to regard the short- 
range part of the two-body interaction, which dominates 
Mq, to alter the statistics only. The long-range part of V2, 
giving the higher moments, modifies the self-consistent 
mean field. Consider for example the electrons in two 
dimensional quantum dots. The two body potential is 
usually taken to be the Coulomb interaction and the con- 
fining potential of the device is modelled by the oscilla- 
tor potential. However, it is expected that the effective 
two-body interaction after averaging over the probability 
densities in the direction perpendicular to the plane will 
be more complicated. Many qualitative features of the 
system may be explained by several choices of the po- 
tential. As in the case of FQHE liquids, we assume that 
the model interaction has a short range part V2sir) and 
a long range part V2i{r). We use the moments expansion 
for the short-range part and neglect the effect of higher 
moments. The self-consistent equation for the density is 
then given by. 



Tn 

P{f) = —p2-{l^-U{r)), r<ro 
ttTi a 



(20) 



= 0, r > ro, 

where the mean TF potential is defined as 

U{r)^VAr)+ J d\p{r')V2i{\f-f'\). (21) 

The equation further simplifies for circularly symmetric 
density. Expanding the potential in partial waves, 

^ oc 

V2i{\r-f'\) = - Vjn{r,r') cosm(6' - 6*'), 

m— 

the TF potential reduces to. 



U{r) = Vi(r) + / r'dr'p{r')vo{r,r'). 



(22) 



In the above equation we have ignored the exchange ef- 
fects which are not important for the long range poten- 
tials. Thus the Eq.(|2^) is the self consistency condition to 
determine the density p{r), and in general is not solvable 
analytically. 

Finally we consider briefly the finite temperature prob- 
lem using the Thomas-Fermi method. We restrict our 
attention to the case where the two body potential is 
extremely short-ranged and regard the system as ideal. 
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The temperature T is expressed in units of the Boltz- 
mann constant, so that it has the dimensions of energy. 
The one-body potential is now temperature-dependent, 
and is given by 



Af 

V{r,T)^V,[r) + ^p{r,T) 

= yi(r)-(l-a)^p(r,T), 



(23) 



where a is the statistics parameter defined by Eq.(15). 
We have assumed that the external potential Vi{r) is 
temperature independent. In the above equation, the 
density p(r, T) for the fermions is obtained from the re- 
lation (including the spin-degeneracy of 2) 



{2'KhY J exp[(p2 /2m* + V - ^l)|T] + l 



(24) 



and the chemical potential is determined by = 
J (Prp(r,T). The p— integration above may be done an- 
alytically, giving 



pir,T) 



m*T 



\n{l + cxp[-{V-p)/T]) 



(25) 



This is inverted to give 

^ = [V + ^P]/T + In (1 - eM-^^'p/^*T) ' 
I m* ^ 



Substituting for V above from Eq.(|2^), we get 



(26) 



In (1 e:^-p{--iih^p/m*T) ' 



(27) 



For a gas in the thermodynamic limit, we set Vi{r) — 
above. Further, the spatial density p may be expressed as 
2/9o, where po is the density for spin- less partcles. Then 
Eq.(27) reduces to the form 



rr. c^—^Po + In (1 - exp{-2nh^po/m*T) ) . (28) 

1 711*1 

Note that this is precisely the equation derived by Wu 
(see his Eq.(23) ) for a two-dimensional gas obeying 
the statistics 



n{e) 



1 



(29) 



w(exp(e — p)/T) + a 
with w{x) satisfying the functional equation 

w"{l + wy-°' ^x^ejip{e~a)/T . (30) 
Here, as in our case, a = 1 corresponds to free fermions. 



We have thus shown that in the large-N limit, ideal 
exclusion statistics may be realised in a system of spin- 
half fermions with very short-range interactions. Note 
that this situation is peculiar to two dimensions since 
both the leading term in the moments expansion and the 
kinetic energy density have the same dependence on the 
spatial density. 
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